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Abstract 

Addition of boundary localised kinetic and Yukawa terms to the action of a 5- 
dimensional Standard Model would non-trivially modify the Kaluza-Klein spectra and 
some of the interactions among the Kaluza-Klein excitations compared to the minimal 
version of this model, in which, these boundary terms are not present. In the minimal 
version of this framework known as Universal Extra Dimensional model, special assump¬ 
tions are made about these unknown, beyond the cut-off contributions to restrict the 
number of unknown parameters of the theory to a minimal. We estimate the contribu¬ 
tion of Kaluza-Klein modes to the branching ratios of —)• pA in the framework 

of non-minimal Universal Extra Dimensional, at one loop level. The results have been 
compared to the experimental data to constrain the parameters of this model. Erom the 
measured decay branching ratio of Bg —)• p~ (depending on the values of boundary 
localised parameters) lower limit on R~^ can be as high as 800 GeV. We have briefly 
reviewed the bounds on nmUED parameter space coming from electroweak precision ob¬ 
servables. The present analysis {Bg —)• p~^p~) has ruled out new regions of parameter 
space in comparison to the analysis of electroweak data. We have revisited the bound on 
R~^ in Universal Extra Dimensional model, which came out to be 454 GeV. This limit on 
R~^ in Universal Extra Dimensional framework is not as competitive as the limits derived 
from the consideration of relic density or Standard Model Higgs boson production and 
decay to W~^W~. Unfortunately, B^ —)• p'^p~ decay branching ratio would not set any 
significant limit on R~^ in a minimal or non-minimal Universal Extra Dimensional model. 


PACS No: 11.10 Kk, 12.60.-i, 14.70.Hp, 14.80.Rt 

I Introduction 


After the discovery of the Higgs boson at the LHC experiment, the new challenge to particle 
physics is to provide a framework in which there exists a natnral Dark Matter (DM) candidate, 
as the Standard Model (SM) itself does not have a snfhciently massive weakly interacting 
particle to be a good candidate for DM. Thns one is compelled to look beyond the SM and in 
this endeavonr, extra dimensional scenarios offer snch a paradigm. Some of the variants of extra 
dimensional theories offer the solntion to the DM Bia pnzzle along with many others like gange 
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coupling unifications ^ and fermion mass hierarchy [1] to name a few. A particular extension 
of the SM needs special attention in this regard. This is known as Universal Extra Dimensional 
(UED) Model where all the SM helds can propagate in 4 +1 dimensional space-time. The extra 
dimension (say, y) is compactihed on a circle (S*^) of radius R [5]. 5-dimensional (5-D) action, 
which has the same held content as the SM, would respect the same SU{3)c x SU{2)l x U{1)y 
gauge symmetry. The 4-dimensional (4-D) effective theory is characterised by Kaluza-Klein 
(KK) towers corresponding to each SM held. The mass of a KK-mode excitation is given 
by = im? + n being the KK-number which is nothing but the discretized momentum in 
^/-direction and m is the 0-mode mass. SM particles have been identihed with the n = 0 mode 
helds in the ehective theory. However there is one caveat, unlike the SM, the 0-mode fermions 
in the ehective theory are not chiral, but are vector like in nature. To get rid of the unwanted 
fermion 0-modes, one needs to impose an extra Z 2 symmetry: y —)• —y on the action. Fields 
which do not have 0-modes are chosen to be odd under this Z 2 symmetry. Rest of the helds, 
having 0-modes, are even under this discrete transformation. The space of y is restricted from 
0 to ttR and is called an /Z 2 orbifold. Two boundary points of the orbifold are also called 
the hxed points as they transform onto themselves under this Z 2 transformation. 

Radiative correction plays mm an important role to cure a highly degenerate KK-mass spec¬ 
trum in this model. These corrections fall into two categories: namely the hnite bulk correc¬ 
tions and more importantly the boundary localised corrections having logarithmic dependence 
on the cut-oh scale A. A very special assumption is being taken in minimal UED (mUED) 
which ensures vanishing radiative corrections at the scale A. In a more general scenario like 
non-minimal UED (nmUED) |8]-|ll], this special assumption is being avoided. Instead, one 
assumes the boundary localised corrections as free parameters of the model. In this article, phe¬ 
nomenology of a particular non-minimal scenario in which kinetic and Yukawa terms involving 
helds are added to the 5-D action, at boundary points, will be investigated. Coefficients of 
boundary localised kinetic terms (BLKT) and boundary localised Yukawa terms (BLYT) along 
with the radius of compactihcation, R can be chosen as free parameters and experimental data 
can be used to constrain them. 

Various phenomenological studies in the framework of nmUED have been made to constrain 
non-minimality parameters from the perspective of electroweak observables [13], S, T and U 
parameters [TTl [T5|, relic density Haim, measurement of decay width of Z-boson to a pair 
of 6-quarks [IB], SM Higgs boson production and decay [IH] and from the LHC experiments 

[201EP. 

Historically any search of new physics beyond the SM has been guided by the effects like preci¬ 
sion electroweak variables like p(T)-parameter, Rf, (Z-boson decay width to a pair of b quarks 
normalised to total hadronic decay width), (forward-backward asymmetry of b quarks at 
Z-pole) etc. Incidentally, all of these electroweak observables are extremely sensitive to the 
quantum corrections. Furthermore, large top quark mass plays a crucial role to amplify these 
quantum effects. In the same spirit, we would like to investigate how a precisely measured 
quantity namely, branching ratio of Bs{d) meson decay to could help us in this endeavour. 

Recent experimental measurement of these branching ratios by CMS [ 22 ] and LHCb [ 23 ] col¬ 
laborations have created some excitements among the high energy physics community, as the 
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Bs —)■ decay has been thought to be one of the harbingers of new physics. Unfortunately, 

the measured value of the branching ratio is more or less consistent with the SM prediction 
[23]. There exists only a one (two)-standard deviation difference between the experimentally 
measured value and the SM estimation of the Bg {B^) branching ratio to /jT pair. However, 
this close agreement of experiment with the SM, can be exploited to constrain the parameters 
of any BSM scenario. In this article we will be performing such an exercise. We will evaluate 
the Br{Bs(d) —t in nmUED framework and compare our results with that from the 

experiment. 

Br{Bs{d) fhe framework of UED has been previously estimated in ref. [25]. How¬ 

ever, presence of boundary localised terms (BLT) in the action would non-trivially change 
masses of KK-excitations and some of the couplings involving KK-excitations in nmUED frame¬ 
work. Thus it would not be a trivial rescaling of an earlier calculation [25] of Bs(^d) To 

the best of our knowledge, this is the hrst estimation of Bs(^d) branching ratios in the framework 
of nmUED. 

In the following section, we will evaluate relevant interactions and vertices in the framework 
of KK-parity-conserving nmUED with a brief introduction of the model. In section 3 and 4 
we will present some calculational details and numerical results respectively. The electroweak 
precision variables have been always very effective in constraining any form of new physics. We 
will briefly review the effect of precision observables (namely the S, T and U) on the nmUED 
parameter space at the end. We conclude in section 5. 


II A lightning review of KK-parity-conserving nmUED 
scenario 


In this section, we discuss the KK-parity conserving nmUED model in brief, with focus on the 
interactions necessary for our calculations. We will need gauge and Yukawa interactions of the 
fermion KK-modes. Furthermore, we will briefly touch the spectrum of Goldstone bosons and 
physical scalars as they play an important role in our analysis as we proceed to present a one 
loop calculation done in Feynman gauge, which is characterised by the presence of unphysical 
degrees of freedom in the particle spectra. 

Let us start with the 5-D fermionic action with BLKT [T^ IT4] : 

S fermion J (fx + rf{6{y) + 6{y - R)]^ D l 

+^riR^Dm-^r + rf{5{y) + 6{y - r] . ( 1 ) 


In 5-dimensions four component Dirac spinors ?/) and ?/) can be written in terms 

of two component spinors [El [13]: 


^dx^y) 


f(t>L{.x,y)\ ((t)^L\x)fdy) 

\XL{x,y)) 



( 2 ) 
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'^R{x,y) 


(3) 



F-matrices satisfy Clifford algebra in 4+1 dimensions: {FmjFjv} = where {M,N = 

0,... 4)0 with the metric gMN = diag(l, —1, —1, —1, —1). The F matrices are dehned as F^ = 
{ 7 ^, ^ 7 ^}, with (/i = 0,..., 3). The 5-D covariant derivative is dehned as Dm = + 

igf^BM, with and gf are the respective 5-D SU{2)l and U{1)y gauge coupling constants. 
^ and y are the corresponding generators. 

Strength of the BLKT for the fermion helds is parametrised by rj which we choose to be 
the same for and T/j to respect the chiral symmetry. We choose equal strengths of the 
boundary terms at the hxed points (y = 0 and y = tiR) which leads to KK-paritj0 conserving 
interactions. We will stick to this convention while dehning the gauge and Yukawa interactions 
in the following. However, in general one can adopt asymmetric boundary terms leading to 
non-conservation of KK-parity. Phenomenology of such KK-parity-non-conserving scenario can 
be found in [HIEI]. 

One can easily obtain the equation of motion governing the dynamics of /l(_r) and gL{R) from the 
above action. With appropriate boundary conditions, solutions are as the following [aiisiin]: 


and 


rn _ n 

Jl — 9r 



cos [mf(n) {y-^)] 


cos 




-sin [mf(n) {y -^)] 


sm 




for n even, 
for n odd. 


9l 


/ 



sin {y- ^)] 


cos 




tR. 


cos {y-^)] 


sm 


'm^(n)TrR 


for n even, 
for n odd. 


(4) 


(5) 


N^, being the normalisation for KK-mode, could be determined form orthonormality con¬ 
ditions: 


lo^dy [1 + rf{6{y) + 6{y - nR)}] f^fl ) dy 1 ^ 

dy [1 + rf{6{y) + 6{y - tiR)}] g^gl J ’ dy gfgl J 

^M,N = 4 corresponds to the j/-direction in space. 

^After restricting the y-direction between 0 and ttR and imposing a discrete Z 2 symmetry on the action 
the residual symmetry of the action is a reflection symmetry along the line y = The transformation of 
KK-excitations under this reflection is characterised by KK-parity. An immediate consequence of KK-parity 
conservation is, in any interaction vertex involving KK-excitations, the sum of KK-numbers must be an even 
integer. 
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and it takes the form: 


N’f = 



1 _L ^ I IL 

4 IT R 


mj^(n) is the n*" KK-mode solution of the following transcendental equations [9l ITT]: 


rfrrifin) 


tan 

cot 


2 


for n even, 
for n odd. 


(7) 


( 8 ) 


We would like to focus on the other relevant interactions (gauge, scalar and Yukawa) which are 
also important for our discussion 




w 

gauge 


(fx 


Sscalar 

Svukawa 


+ ry{6{y) + 5{y - 7rR)}W^^^W;, 
d^x [ {D^^y (Dm^) + r^{S{y) + 5{y - nR)} (D^d>) 

d^x\x^ + Ty {5{y) + 5{y - 'KR)}\l(j)L^XR + h.c. 


(9) 

( 10 ) 

( 11 ) 


WIj^ = {OmWI^ — BnWIj — represents the SU{2)l {i = 1,2,3) field strength 

tensor. $ is the Higgs doublet and <h satisfy the condition $ = with being the Pauli 

matrix, ry and are the coefficient of BLKT for gauge and scalar fields respectively while 
Ty represents the coefficient of boundary terms for Yukawa interactions. denote the Yukawa 
interactions strength for the third generations in 5-D theory. Further, the 5-D gauge coupling 
is related to the 4-D coupling g 2 through. 


gl = g 2 




Appropriate KK-expansion of the fields which are involved in the above actions can be schemat¬ 
ically written as: 

W^(x,y) = ^lvy)(i)o”(„) ; Wl(x,y) = ^ in‘'”A)6”(!/). (12) 

n n 

and 

^x,y) = J2^^-\x)h-{y). (13) 


Note, that ?/-dependent profiles for the and must be different in the view of the fact 
that KK-tower must have a 0-mode, while the later shouldn’t have such a mode in its KK- 
tower. In general, due to electroweak symmetry breaking, the eigenvalue equation for the gauge 
boson contains a term proportional to (r<^ — ry)[26]. Consequently, KK-solutions of the gauge 
bosons are in general, different form Eq. IH To avoid unnecessary complication we set = ry, 
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which we will keep in the rest of our analysis. Consequently, gauge boson KK-excitations have 
masses my(n) (= which satisfy the same transcendental equations given in Eq. |8]with r/ 

replaced by ry (= r^). Furthermore, this assumption helps us to fix the gauge properly in this 
non-trivial scenario. For a detail discussion of the and ^/-dependent profiles we refer to 
the ref. [2S]- 

One must supplement the action for fermions and gauge bosons with the gauge-fixing action 


qW _ J- 

“ 'T 
sy 

(14) 

Here, is the 5*^ component of charged SU{2) gauge boson, while (j)'^ is the ^ 

component of the Higgs doublet. Mw is the 0-mode IF-boson mass and is connected with 
physical gauge fixing parameter ^ (with values 0 (Landau gauge), 1 (Feynman gauge) or oo 
(Unitary gauge)) via, 

^ = 4{1+ ry (5(|/)+ 5(|/-7ri?))}. (15) 

Using Eqs. [9l [TOl and Eq. [THwe can write the bi-linear terms involving the KK-modes of 
and in gauge as: 


d^x 


+ rv {5{y) + 5{y - vri?))}) 


C 


y/4(n)±^(n)± 




0(0-) 


f 

0i(l - Cl^wruyin) 


"iyw + / V ) ’ 


(16) 


In addison to Goldstone bosons (with mode mass square ^(M^ -|- m^(„))); 

Myr(n) 

we have additional physical charged Higgs pair (with mode mass square (M^ -|- m^(„))): 


Mw(n) 


The helds H^^("')=*= and have the same mass eigenvalue Myr(n) = ^and 

’t-Hooft Feynman gauge (.^ = 1), G'("’)=*= also correspond to the same mass eigenvalue. 


m 


Let us now examine the mixing of the quark sector. This mixing is only important for top 
quarks as it is driven by the Yukawa coupling. Using the modal expansions for fermions given 
in Eqs. [2] and |3] and substituting these in the actions given in Eq. [1] and Eq. [11] we can End the 
bi-linear terms involving the doublet and singlet states of the quarks. The mass matrix for 
KK-level is as the following: 
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[(pL ) y -rrifin) 



(m) ) + h-C., 

xk ' 


(17) 


where m* is the top quark mass and mj(n) are the solutions of transcendental equations given 
in Eq. |H1 The overlap integrals and ^ 2 ^) are of the form: 

/f I IL \ r-KR 

A”” = () X / dy[i + r„{S(y) + Hv - sS/?, 


and 


nR 


' 1 _|_ \ /"Tri? 

^nm _ / -kR 


1 -L 
^ TTiJ 


X dy g^fl. 


For both the cases of n = m and n ^ m the integral is non zero. But for Vy = Vf, 

this integral equal to 1 (when n = m) or 0 (n 7 ^ m). And the integral ^ 2 ^ is non zero only 
when n = m and equal to 1 in the limit Vy = Vf. Using this equality (fermion and Yukawa 
BLT) condition we can easily avoid the mode mixing and construct simpler form of the fermion 
mixing matrix. In the rest of our analysis we will stick to the choice of equal Vy and rj. 

The resulting matrix (given in Eq. [T7)l can be diagonalised by following bi-unitary transforma¬ 
tions for the left- and right-handed helds respectively: 


_ / COS (y.tn Sm O^tri 

^ '-sinatn cos atn 


_ I COS cutfi sm C^tn 

^ I sin atn - cos atn 


(18) 


where atn[= |tan ^ mixing angle. The gauge eigen states '^L{x,y) 

and \k/j(a:,|/) and mass eigen states and are related by the following relations: 


y ( 77 ) rTil(77.) • rT-i2(77-) 

(Pl = cos atnT^l - sm atnT^l % 
= cos atnT^p^ + sin atnT^p\ 


= sin atnTljp^ + cos atnTp''\ 

X^R = sin atnTlP^ - cos atnT^p ^. (19) 


The eigen states are degenerate in mass, values of which are given by: rUrpHn) = m.^ 2 (n) = 

Feynman rules (in ’t-Hooft Feynman gauge) which follow from the above action and necessary 
for our calculation are listed in Appendix C. One can see that some of the interactions listed 
in Appendix C have been modihed with respect to their UED counterparts by some overlap 
integrals. These overlap integrals are generated while arriving at the 4-D effective action from 
5-D one by inserting the appropriate ^/-dependent prohles in the (5-D) action and integrating 
over the ^/-direction. Here we present the overlap integrals which will frequently occur in our 
calculation. 
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Interaction of a O-mode fermion with a KK-fermion mode) and a KK-gauge boson 

(jrith jYiode): 




ry \ 


X 



dy [1 + rf{6{y) + S{y - vri?)}] , 

.ii ( 20 ) 

dy [1 + rf{5{y) + 5{y - nR)}] . 


Interaction of a 0-mode fermion with a KK-fermion mode) and the KK-mode of 
^th component of a gauge boson: 

pTrR 

/ dy b^glfl , 


jfnm 


= \ nR\l + 


ry \ 
-kr) 


X 




dy b’"figl. 


( 21 ) 


In Eq. [20] a”^ is the wave function KK-mode) for gauge held, with Lorentz index 
y. While in Eq. |2T] 6™ is the wave function KK-mode) for 5*^ component of the 
corresponding gauge held. 


• Yukawa interactions of fermions and scalars: 

dy [1 + r,y{y) + 5(y - irB)}] , 

dy |1 + ri{6(y) + d{y - irB)}] . 

In the above, h’^ is the wave function (ml" KK-mode) for scalar field. 



( 22 ) 


Now for, = ry we have and for Vf = Vy, ly'^ = for n = m let us call it /{^, and 

is given by: 


/r = 2 i 


'1 + 2 V 
^ ttR 


1 + 


ttR 


Tfra, N 

1 _L ^ I 11. 

4 TTiJ 


1 + 


^v'^v(n) 


I rv_ 
ttR 


m 


y(«) (rf - ry) 




7tR 


(23) 


and we denote ly"^ = Let us call it for n = m and is given by: 


/l+^ 

1 


1 

lTly(n) TR j{‘n) 

1 

1 


4/1 ++ -f 

^/i 1 1 r-vv 

LV ”'" 4 “’"ttR-I 


tiR 


( 24 ) 

































Rf Rf 

Figure 1: Variation of the overlap integrals II (left panel) and I\ (right panel) dehned in Eqs. |23] 
and [21] with Rf for six different values of Ry- 


These two overlap integrals play a crucial role in our analysis. Many of the couplings used in 
the following calculation are hallmarked by the presence these integrals which can be explicitly 
seen in the Feynman rules listed in Appendix C. First of all, in the limit rj = rv, both of these 
overlap integrals become unity, setting all the nmlJED couplings equal to their UED values. For 
the purpose of illustration, variations of the integrals I\ and {n = m = 1) with rj/R{= Rf) 
and ry/R(= Ry) have been presented in Fig. [TJ Left (right) panel shows the variation of II 
(J 2 ) with Rf for six different values of Ry. We have presented these integrals for both positive 
as well as negative values of Ry and Rf. Let us hrst discuss the pattern of these integrals for 
positive values of Ry. In this case, while l\ increases with Rf,l\ initially increases before it 
starts decreasing with Rf. On the other hand if Ry is negative, /| is almost independent of 
Rf, but Jg steadily decreases with Rf over the entire range of its values used in Fig. [H At 
this point we would like to make some comments which could help the reader to understand 
the nature of the plots. The hrst two factors (one under the square root sign and the other in 
the hrst square bracket involving rf only) decreases with i?/ in both of these integrals. Among 
the remaining factors, term proportional to (in l\) slowly diminish in magnitude with 

Rf. The remaining factor proportional to {Rf — Rv) increases with Rf. When Ry is negative 
(= —1), due to an interplay between the term proportional to m^(i) and remaining factors, II 
remains almost independent with Rf. This explains the qualitative departure of the overlap 
integral II for negative and positive values of Ry. Furthermore, as rrifii) is a decreasing function 
of Rf, in II the presence of a factor proportional to rnymuify) (instead of in Jj^) ohers 

additional damping to the magnitude of Jg. And apart from some lower values of i?/, J 2 always 
decreases with Rf. 

Before going into the main calculation, let us comment on the range of values of BLT parameters 
used in our analysis. In general BLT parameters can positive or negative. A careful look into 
the Eq.[7| would reveal that, for Vf/R = —tt the 0-mode solution becomes divergent. And 
beyond rf/R= —vr the 0-mode helds become ghost-like. Any other values of BLT parameters 
greater than —vr are acceptable. However, as BLT parameters change from positive to negative 
domain, corresponding KK-masses would increase thus diminishing the values of loop functions 
described below. This in turn decrease the magnitudes of the branching ratios of our interest. 
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Ill in nmUED 

Stage has now been set to discuss the subject of our attention, namely the branching ratio of 
Bs(d) into pair. The effective Hamiltonian for i?s(d) —t decay is given by: 

Cjr 

‘H.s = -^^^^^^V,lV„(_iyriyY(x,,rf, rv,R~') [fr7„ (1 - 75)s(<i)] 1^7'“ (1 “ Ts)^] + h.c. , 

(25) 

where Gp is the Fermi constant, Vij are the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele¬ 
ments, a is the hne structure constant, r]Y is the QCD factor and is the Weinberg angle. 
The function Y{xt,rf,rv,R~^) is the total contributions coming from Z-penguins, self-energy 
and box diagrams: 


Y{xt,rf,rv,R ^) = C{xt,rf,rv, R ^) + B{xt,rf,rv, R ^), 


(26) 


where C{xt,rf,rv, R originates from penguins and self energy diagrams (Fig. |2] and Fig. [3] 
respectively) and is dehned as: 


OO 

C{Xt,rf,rv,R~^) = Co{Xt) + '^Cn{Xt(n),Xu(n))- (27) 

n=l 


Here, Co{xt) is the SM contribution: 


Coixt) 


8 


xt- 6 3xt + 2 
Xt-1^ {xt- ly 


\nxt , 


(28) 


while the second term represents the total KK-contribution that is being calculated from the 
Z-penguin diagrams and self-energy diagrams (Figs. [2] and [3]) originating from the nmUED 
framework. The function Cn{xt(n),Xu{n)) is dehned as: 


Gn{Xt(ji'^i X^^^y 7^(^ti(n))5 (29) 

where F{xt(n)) and F{xu{n)) arises respectively, from the contributions of the and 

Tu^^\ modes. The function F{xt(n)) takes the form: 


8 

7^(^t(n)) ^ ^ Rj iXt(n )) T 

i=l 



- sin^ ) 5]^ /\Si{Xt(n))- 


2 = 1 


(30) 


Moreover Fj denote the contributions those are coming from Z-penguin diagrams (1-8 in Fig. [2]) 
and A5j represent the contributions of self-energy diagrams (la-IIb in Fig. [3]) which are nec¬ 
essary for calculation of the electroweak counter terms. While the function F{xu(n)) is given 
by: 


7^(^ii(n)) 7^(^t(n)) 


(31) 


xt—^0 
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F{xu(n)) are subtracted in the above (Eq. [2nD to take into account of the contributions of KK- 
excitation of first two generations of quarks in Z-penguins and self-energy diagrams, exploiting 
the GIM-mechanism. F{xt{n)) in nmlJED framework is given by the following expression which 
is drastically different from that of UED due to the presence of integrals /{*and : 


F{Xt{n)) = 


- InM^t^) -- + hq {xt^n)) \ - -:Xt{n)hq {xt(n)) 


m 


fir.) 




d" 


n\2 


1 1 fin) 

^^ wi '^) ~ 5 ~ ^t{n)hw (^i(n)) (-^1 ) ~ 2 ^^ (^^('^)) j ^2 ) 


Wi^) 


+ 


1 


1 


\Q I r) F hq [x^n)) f {^2 ) QXt(n)hq (xt(n)) 


1 


m 


fM 


m 2 


-f 


4 

mt 




^t(n) M^ 


(Fi 


n\2 




-f- 


16 

1 


1 f l + x,(.) ^ 2 x 2 („)lnx,( 


2)1 Xi(ji^ (1 Xf(jf) 


n/(") 


{II 


n\2 


1 






f Xt{n) (1 ®t(n)) 


2 f 


(32) 


Expressions for hq and are listed in Appendix A (Eqs. A-9 and A-10). 



Figure 2: Z-penguin diagrams contributing to the decay of -Bs(d). 
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Figure 3: Self-energy diagrams contributing to the decay of Bs(d)- 


B{xt,rf,rv, R in Eq. [2n]is given by: 

CO 

B{xt,rf,rv,R~^) = -Bo{xt) + '^Bn{xt{n),Xu(n),x^^ri))- (33) 

n=l 

Bq^Xi) again represents the SM contribution to the box diagrams: 


Bo{xt) 


1 

4 


Xt Xt 

1-xt {xt - ly 


\nxt , 


(34) 


and the second term denotes the total KK-contribution that is evaluated from the box diagrams 
(Fig.i]). Bn{xt(n), Xu(n), IS deEued as: 

Bn{Xt{n)-i Xu(n)-i Xi,(jiy H {Xi(^yi)i H (^Xu(^ji) j Xi/(^ji)') ■ (35) 

In nmlJED framework, expression for H{xt{n)iXy(n-)) is given as the following: 


H{Xt{n),X^(^n)) = {Ilf 

1 


1 ~ 

7Tp R{Xt{n)-!Xu(n)) + X R{Xt{n)-:Xv{n)) 


1 


+ 


0 2 \ m ^ 2 


2 


m 






8 ““«») 


m 


y(^) 


-mj U{Xt(^n),Xuin)) 


1 

^ <(.) 


1 


U{Xt(n)^Xy(^n)) IQ 


4 +mimy^)\ U{Xt(n),Xu{n)) 

WM \ / 


(36) 


and H{xu(n),Xu{n)) is dehned through: 


H {Xu{n)) Xy(Yl)) H {xt^n) 1 Xiy^^jif 


(37) 


xt^O 


The expressions for the functions U and U can be found in Appendix B (Eqs. B-7 and B- 8 ). 
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The definition of Bn (given in Eq. 1531) are again in order to take into acconnt the contribntions 
from KK-excitations of first two generations of qnarks into the calcnlation. 




Tin 




Uj , Cj 






Figure 4: Box diagrams contributing to the decay of i?s(d). 


Note that Xt = 


'W. 


■, and Xt(n) = , rUfin) and my(„) are the solutions of Eq. [8] 


corresponding to the BLKT parameters Vf and ry respectively. 


Finally one can express the branching ratio of i?s(d) —)■ in terms of Y and other parameters 

as: 


Br{B, 


's(d)^n^^ ) - nBs{d))^(- -. 2n 1 

TT ydvrsm 0 ^/ 


Fr rrlmB 

t>s{d) A* 


(d) 


1 - 4- 


m. 


m 


nt,rf,rv,R-l, 


B. 


(d) 


(38) 


where, is the i?s(d) meson decay constant and is the life time, and are 

the masses of muon and i?s(d) meson respectively. 

We would like to point out that the total contribution from UED is obtained by summing CnS 
and BnS over KK-modes (starting from n = 1). It can be checked easily that these quantities 
tend to vanish when R~^ oo, showing the decoupling nature of the new physics in this 
case. While estimating the contribution from the KK-excitation, in the one loop mediated 
amplitudes, we have used only those interactions which couples a 0-mode field to pair of KK- 
excitations with same KK-number. However, in a manifestly KK-parity conserving theory, 
there can be non zero interactions among a 0-mode field with two different KK-excitations, 
say, with unequal KK-number n and m, modulo n -|- m is an even integer. We have explicitly 
checked that 0 — n — m interactions are suppressed with respect to the 0 — n — n interactions 
[18] and we neglect those while estimating the Bs{d) interactions. 
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In our calculation, we have neglected the masses of the external quarks and while evaluating 
the Z-penguins and box diagrams, momenta of the external legs have been assumed to be 
vanishing. 


IV Results 


We are now equipped to present our results. Before we proceed further let us present the 
experimentally measured value m of Br{Bs{d) —t ) along with the SM prediction [Ml ITT] 

for the same in Table [1] Central values of experimentally measured branching ratios of Bg {Bd) 
are nearly one (two)-standard deviation away from its SM prediction. 


Branching ratio 

SM prediction 

Experimental value 

Bg —)■ 

(3.66 ±0.23) X 10“® 

{2.S1Z) X 10-» 

Bd 

(1.06 ±0.09) X 10“^® 

(3.9li;^) X 10“^® 


Table 1: SM prediction for Bs(^d) branching ratios and their experimental measured 

values. 


We will now discuss the method, that has been used in the following, for the numerical estima¬ 
tion of Br{Bs —)■ We will follow similar procedure for the estimation of Br{Bd —)■ 

Numerical value of the Bg branching ratio, quoted in Table H] includes higher order corrections 
(NNLO-QCD and NLO-EW corrections [23]) to the branching ratio in the framework of the 
SM. The total decay amplitude in the UED framework is a sum of contributions from different 
KK-levels (including the contribution from 0-level which is nothing but the SM). We expect 
that higher order corrections to the decay amplitude originating from higher KK-modes would 
be in the same ballpark that of the SM (0-mode). So to (approximately) take into account the 
effects of higher order corrections into our result we have estimated Br{Bs —)■ /r^/i“) in UED 

asj^ 

Br{Bg —)■ )uED = Br{Bg —)■ )sm 

Here, Y total contribution (originating from 0 and higher KK-modes) dehned in Eq. l26l while 
Yq{= Cq + Bq) is only the 0-mode (SM) contribution. Expressions for Cq and Bq are given in 
Eq. |28]and Eq. [Ml respectively. For numerical estimation of the above branching ratio, we have 
used the central value (3.66 x 10“®) for Br{Bg —)■ following the ref. [27] . 

In the above, we have expressed Br^Bg fJ^~)uED in terms of Bri^Bg —)■ Con¬ 

sequently, the relevant input parameters for numerical evaluation of the above branching ratio 

^This method is somehow similar to find a A:-factor (the ratio of branching ratios with higher order correction 
and without any correction in our case) in the SM and then multiply the leading order branching ratio in UED 
by this fc-factor. In absence of a full higher order correction in UED framework, this approximation is not very 
bad, as the higher KK-modes of t-quark and lU-bosons, running in the triangle or box diagrams, have similar 
couplings to the SM particles in the initial and the final state. 
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are top-quark and ly-boson masses respectively, which are necessary for evaluation of the func¬ 
tions Yq and Y. We use the values of these parameters following the ref. [28]. Values of these 
parameters used in ref. [ 23 ] are same as the values in ref. [280. 

Any contribution to the i?s(d) —)■ branching ratio from the SM and beyond the SM goes 

into the function Y as defined in Eq. [381 The dominant contribution to Y comes from the 
diagrams (3), (4), (6) and (7) of Fig. [2] (Z-penguins); (Ib) and (Ilb) of Fig. [3] (self-energies). 
These amplitudes are proportional to g 2 ^t^ where Xt being the top-Yukawa coupling. Box 
diagram contributions are not significant in comparison to the penguins [§. One must perform 
a sum over the KK-modes while evaluating the C and i?-functions, which contribute to Y. 
In view of a recent analysis correlating the SM Higgs mass and cut-off scale of UED model 
[29], we restrict ourselves upto 5 KK-levels. Previous practice was to use 20-30 KK-levels 
while summing up the contributions from KK-modes. However, we have explicitly checked 
that numerical results would not change to a signihcant level as the sum over the KK-modes, in 
this case is converging. At the end of following sub-section, we will present a table comparing 
the Bs branching ratios calculated with 5 KK-levels vis-a-vis 20 KK-levels in support of our 
assumption. 

The main results of our analysis i.e. variation of the branching ratio {Bg —)■ with R~^ 

for several values of Rf are presented in Fig. |6l Same have been presented for B^ in Fig. [71 
Six different panels show the dependence of branching ratios with Ry- We will be discussing 
the implications of these results in the framework of nmlJED. But before delving into that, we 
would like to discuss the consequences of the result in the framework of UED itself. 

IV. 1 Reviewing the limits on R~^ in UED framework 

To start with we will focus on the Bg^d) meson branching ratio to gY in the UED frame work. 
The expression for F” and defined in Eqs. [32] and [36] would result into their UED forms once 
we set ry = Vf = 0. In this limit the overlap integrals and become unity and the masses 
of the KK-excitations in KK-mode become equal to nR~^. We check that our expressions 
of C and B in this limit agree exactly with the one in ref. [250 . In Fig. [5] we have plotted 
the Hs(d) —t gY branching ratio with R~^ in UED framework (black dashed curve in each 
panel). Steady decrement of the branching ratio with increasing R~^ is an artefact of increasing 
KK-masses with R~^. This line intersects the horizontal line (light shaded) (4.2 x 10“® for Bg 
and 7.1 X for Bd) which correspond to the 95% C.L. upper limit on the experimentally 

measured branching ratio of Bg (left panel of Fig. [5]) and Bd (right panel of Fig. [5]). Thus the 

'^We have used Mw = 80.38 GeV and mt = 173.21 GeV. 

®As for example, in nmUED model, for R~^ = 1 TeV, Rf = Ry = 1 after summing over 5 KK-levels 
C = 0.9116 while B = 0.1856. For R-^ = 2 TeV, Rf = Ry = 1, C = 0.8915 while B = 0.1846. The smallness 
of the box diagram contribution in comparison to the penguins, has also been pointed out in ref. [25] 

®It has been shown in [30], the KK-sum in case of one loop calculation (like the one in the present case) in 
UED with one extra space like dimension is always converging. 

^The authors of ref. [25] have not considered any radiative corrections to the KK-masses, consequently mass 
of a KK-excitation is nR~^ in their analysis. 
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0.2 0.4 0.6 0.8 1.0 

(TeV) 


Figure 5: Left (right) panel: Variation of the branching ratio {Bs{Bd) —)■ with R~^ in 

UED model. The black dashed line in each panel represents the UED curve (corresponds to 
ry/R = Rv = 0 and rf/R = Rf = 0 case). The dark shaded horizontal line represents the SM 
prediction for the above branching ratio, while the light shaded horizontal line represents the 
95 % C.L. upper limit on the experimentally measured Bs{d) branching ratio to muon pair. 


R~^ values corresponding to these intersections would give rise to a 95% C.L. lower limit on the 
value of R~^. One can see this limit from Bg decay comes out to be 454 GeV. While deriving 
this limit, we have neglected the theoretical error on Br{Bs —>■ One can easily see from 

Table [H that theoretical error on the SM estimate of this branching ratio is around 6% and we 
expect that theory error in the UED framework will be of the same order. And this is small 
compared to the experimental error (nearly 25 %) on the measured value of the branching 
ratio. Our results presented in Fig. El thus correspond to a vanishing theoretical error on the 
Bg branching ratio. 

The lower bound on R~^ which results from B^ decay {R~^ ~ 60 GeV) is not so promising. To 
put our observation into the proper context we would briefly mention the bounds on R~^ that 
result from other processes. Consideration of {g — 2)^ [31], p-parameter [32], FCNC process 
[25l [33l |3l] and electroweak observables like R^ [HI |35] would result into a lower limit on R~^ 
which is in the ballpark of 300 GeV. Consideration of radiative 5-meson decay puts a lower 
limit of 600 GeV on R~^ [36]. While the projected lower limit on R~^ using tri-lepton signal 
at 8 TeV LHC is 1.2 TeV [37] . 


IV.2 Bounds from the nmUED 

We are now in a position to discuss our results in the context of nmUED. Numerical results 
from the calculation presented in the previous section are presented in Figs. |6] and [7] for Bg 
and Bd meson branching ratios to muon pair respectively. In each of the hgures six different 
panels represents the results for six different choices for gauge BLKT, Ry. While in each panel 
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different line represents the variation of the branching ratio with R ^ with different valnes of 
fermion BLKT coefficients, Rf. 

The monotonic decreasing natnre of the branching ratio with increasing R~^ is the same as in 
UED and has been explained in the above. The dependence on the other parameters conld 
be easily understood as the following. For a given masses of KK-excitations of fermions, 
gauge bosons and scalars decrease with increasing Rf and Ry- This in turn, increase the decay 
width of Bg and to . Apart from KK-masses the other determining factors in our 
calculation are the overlap integrals dehned in Eqs. 122] and |2H We have already noted their 
dependencies on Rf and Ry. 

One can see that Bg —)■ branching ratio increases with the increment of both of these 

BLKT coefficients. Dependence on Ry is mild while it is more sensitive to any change of 
Rf. This can be understood again by looking at the interactions used in this calculation 
listed in Appendix C. A careful look at the Feynman rules reveals that /{* and modify the 
interactions of 3^'^ generations of quarks with charged-Higgs scalar and 5*^ component of IK- 
boson respectively. The hrst one being top-Yukawa, is dominant over the second one, which is 
SU{2) gauge interaction. Thus which appears with top-Yukawa, has a better control on the 
Bg —)■ /r■*“/!“ amplitude, which increases both with Rf and Ry. Dependence of Br{Bd —)■ 
on R~^ and BLKT parameters can be explained in the same way. The decay of Bd to R^R~ is 

suppressed with respect to the Bg by a factor of (yyf^ ■ Finally we would like to add that for 
negative values of BLKT parameters, above branching ratio is always less than the value for the 
same in the UED model. This could be easily accounted by the suppression of decay amplitude 
due to heavier (than in the UED case) KK-excitation in nmUED framework for negative values 
of BLKT parameters. 

Finally in Fig. [H] (right panel), the exclusion region in Rf — R~^ plane, has been presented for six 
different choices of Ry. The region under a particular line (corresponding to a hxed value of Ry) 
has been excluded at 95% C.L. by comparing the experimentally measured Br{Bg —)■ /i^/i“) to 
its theoretical prediction in the framework of nmUED. This plot is a summary of the results 
presented in Fig. [6l The lines represent the contours of constant branching ratios of Bg —)■ R^ 
corresponding to the 95% C.L. upper limit (4.2 x 10“®) of its experimentally determined value. 
Features of this contours can be easily understood in the light of our discussion above. Higher 
values of R~^ would increase the KK-masses hence pulling down the decay width (and branching 
ratio). To compensate this one must increase Rf and Ry, which control the decay dynamics 
in two ways. First of all, these would pull down the masses and Rf would also increase the 
couplings through the overlap integral IR While an increasing Ry would decrease but 
increase IR Overall, with increasing Ry, decay width mildly increase in contrast to Rf, which 
would increase the decay width more sharply with its increment. 

We would like to show the difference in the derived lower limits on R~^ due to taking into 
account 5 KK-levels vis-a-vis 20 KK-level while performing the KK-level summation. One can 
see from Table [21 that the lower limits on R~^ (in TeV) derived from our analysis is not very 
much sensitive to the number of KK-level used in the summation. 
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Figure 6: Variation of the branching ratio {Bg —)■ with R~^ for several values of Rf = 

Tf/R. The six panels correspond to different Ry = ry/R. The dark shaded horizontal line 
represents the SM prediction for the above branching ratio, while the light shaded horizontal 
line represents the 95 % C.L. upper limit on the experimentally measured Bg branching ratio 
to muon pair. 
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Figure 7: Variation of the branching ratio (Bd —>■ with R~^ for several values of Rf = 

Tf/R. The six panels correspond to different Ry = ry/R. The dark shaded horizontal line 
represents the SM prediction for the above branching ratio, while the light shaded horizontal 
line represents the 95 % C.L. upper limit on the experimentally measured Bd branching ratio 
to muon pair. 
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Rv 

= -1 

Rv 

= 0 

Rv 

= 1 

Rv 

= 2 

Rv 

= 3 

Rv 

= 4 

Ri 

5 KK-lcvcl 

20 KK-lcvcl 

5 KK-lcvcl 

20 KK-level 

5 KK-lcvcl 

20 KK-level 

5 KK-lcvcl 

20 KK-lcvcl 

5 KK-level 

20 KK-lcvcl 

5 KK-lcvcl 

20 KK-lcvcl 

-2 

0.333 

0.357 

0.246 

0.260 

0.119 

0.122 

0.077 

0.081 

0.045 

0.050 

0.035 

0.039 

-1 

0.361 

0.378 

0.280 

0.300 

0.191 

0.192 

0.146 

0.150 

0.119 

0.120 

0.101 

0.105 

0 

0.368 

0.380 

0.454 

0.460 

0.412 

0.418 

0.392 

0.402 

0.383 

0.392 

0.373 

0.384 

1 

0.381 

0.399 

0.553 

0.567 

0.586 

0.603 

0.588 

0.610 

0.590 

0.617 

0.593 

0.621 

2 

0.399 

0.406 

0.617 

0.624 

0.672 

0.686 

0.682 

0.700 

0.693 

0.717 

0.698 

0.723 

3 

0.403 

0.409 

0.661 

0.664 

0.734 

0.745 

0.752 

0.765 

0.766 

0.786 

0.770 

0.792 

4 

0.410 

0.416 

0.694 

0.701 

0.780 

0.783 

0.800 

0.810 

0.820 

0.834 

0.825 

0.838 


Table 2: Lower limits on R~^ (in TeV) derived from Bg decay branching ratio to for 

different values of input parameters showing the insensitivity on the number of KK-levels in 
summation (see Eg. [271 and Eg. [551). 


IV.3 Electroweak precision constraints 


Before we conclude, it would be relevant to discuss the constraints on the parameters coming 
from our analysis vis-a-vis electroweak data, which remains very instrumental in constraining 
any new-physics beyond the SM. Electroweak precision constraints on nmlJED parameters have 
been discussed previously in [38l [39] , while authors in ref. [13] restricts nmlJED parameters from 
the consideration of Z-mass. However, in both of these cases nmUED action used for analysis 
are slightly different from that of ours. Ref. [38] , while presenting their electroweak results used 
equal BLKT parameters for all the held along with a bulk mass term for the fermions. While 
the authors in the ref. [13] , have considered several choices for the BLKT parameters. The case 
closest to our approach is where they have used egual BLKT parameters for Higgs and gauge 
bosons and set the BLKT parameter for the fermions eguals to zero. Conseguently, it would 
not be meaningful to directly apply the constraints derived in the above analysis to our case. 
So we have redone the analysis following the approaches of refs. [131 EH] , however applied to 
our case. 


In nmlJED model, corrections to Beskin-Takeuchi parameters S, T and U as well as to Z 
mass appear through the correction to the Fermi constant, Gp at tree level, which is in stark 
contrast to the minimal version of the UED model where these correction appears through one 
loop processes. These guantities are modihed by the correction of Fermi constant Gp which is 
determined from muon decay, i.e. a four-fermion process. The corrected Fermi constant Gp 
can be decomposed into two parts as: 


Gp = 


+ SGf, 


(40) 


with G^p is simply the contribution from 0-mode exchange, while 5Gp denotes the sum of 
the contributions from all non zero (even) KK-modes. The effective Fermi constant can 
be expressed as the following. 


^0 _ 92 


n>2 




1 


(41) 
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where, 


Inoo = + [l + rj{5(|/)+ 5(|/-7ri?)}]aVL/E, (42) 

^ a /(1 + ^) (rf-ry) 

(1 + ^) ./TTlS^lT^ 

V “'" 4 “■“ ttK 

A similar expression for SGp has been presented in ref. [I3]. The above expression for 6Gf, 
agrees with the same in ref. [13] in the limit —)• 0. 

Note that the above integral becomes zero when Vf = ry, i.e., there is no correction of Fermi 
constant for this eqnality condition. 

nmlJED contributions to the S, T and U parameters can be written, following the approach of 
ref. [38l |39], as: 




nmUED 


= 0, T, 


nmUED 


16Gi 
ex G p 




4sin^ 5Gf 


nmUED 


a 


Gi 


(43) 


One can now compare the predictions from nmUED model with the experimental values given 
in the ref. jlO] 


,Snp = 0.05 ± 0.11, Tnp = 0.09 ± 0.13, Unp = 0.01 ± 0.11, (44) 

with input Higg mass mh = 125 GeV and top quark mass nit = 173 GeV. This would constrain 
the parameter space of nmUED model. 

Before presenting the results of the above analysis, we would briefly discuss how 5Gf modihes 
the Z-boson mass at the tree level [T3|. This would be evident once we express the mass of Z- 

boson in terms of a, Mw and Gf '■ Mz = _ rppg corrections to Z mass at tree level, 

in the framework of nmUED, would creep in through the corrected Fermi-constant, Gf- One 
can compare this with the so called tree level Z-mass defined by the relation 
following the ref. [TB]. Here, 5^^~^°°^'^mz stands for the 1-loop correction to the Z-boson mass 
in the SM. 

We present our results in Fig. |H] (left panel) in terms of 95% G.L. lower limit on R~^ for several 
values of BLKT parameters. Region of the parameter space below a particular line has been 
ruled out from the consideration of T and U parameters. Furthermore, we have checked that 
the tree level value of Z-boson mass in the framework of nmUED, resides within the 95 % G.L. 
allowed tree level value [13] for the entire range of parameter space that we have considered in 
this work. 

To compare the limits derived from our analysis {Bg —)■ with that from electroweak 

precision test, we present the allowed range of parameter space both from Bg decay and elec¬ 
troweak precision data, in Rf — R~^ plane for six different choices of Ry in Fig. |8] (left panel). 
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In the right panel of Fig. [HI we present the allowed parameter space, on an expanded scale, 
only from the Bg decay analysis for sake of clarity. Each curve (in the left and right panel) 
corresponds to a particular value of Ry, each point of which represents a lower limit of R~^ in 
TeV corresponding to a particular value of Rf. So that the area under a particular curve has 
been dis-allowed by either from precision constraints or from Br{Bs —)• It is evident 

from the plots that for a given value of Ry, higher values of R~^ are being ruled out from 
electroweak data for lower values of Rf, while Bg branching ratio would do better than the 
electroweak data in excluding values of R~^ for higher values of Rf. Finally, we mention while 
passing that the all of the parameter space (whether allowed or disallowed from electroweak 
data and Bg decay) is consistent with the Z-boson mass following the ref. [T3] . 




Rj Rj 

Figure 8: 95% C.L. exclusion contours in Rf — R~^ plane for six different choices of Ry from 
branching ratio of Bg —)• decay and from constraints of T and U parameter. The area 

below a particular curve (hxed Ry) has been excluded at 95% C.L. 


It has been evident from our analysis that the lower bounds on R~^ derived with negative values 
of BLKT parameters are not very interesting and the above values of R~^ have been already 
excluded from the consideration of electroweak data, which can be seen from the left panel 
of Fig. [HI The weak nature of the bounds on R~^ derived from Bg —)■ branching ratio 

for negative values of BLT parameters, could be accounted by the higher masses of the KK- 
excitations, diminishing the contributions of penguins and boxes to the total decay amplitude. 
Furthermore, we have restricted our choice of BLT parameters upto 4. A careful look at the 
left or right panel of Fig. [H] (or Table [2]), tells us that lower limits on R~^ are weakly sensitive 
to the values of Ry. So we expect that higher values of Ry would keep the lower limit on R~^ 
in the same ballpark. 
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V Conclusion 


We have calculated the contribution of KK-excitations in the framework of non-minimal Uni¬ 
versal Extra Dimensional model to the branching ratio of Bs{d) at one loop level. 

Non-minimal UED is hallmarked by the presence of boundary localised kinetic and Yukawa 
terms along with an SM like action and held content however dehned in a 4 -|- 1 dimensional 
bulk. Boundary localised terms parametrise the unknown radiative corrections to the masses 
and couplings in the full 5-D theory. Presence of boundary terms modify the couplings and 
mass spectrum of KK-modes in the 4-D effective theory in a non-trivial manner. 

To put our discussion into a context, we must remind that in UED, masses of the KK- 
modes are nR~^ {R~^ being the compactihcation scale). Interactions among the different 
KK-excitations are very similar to their SM (0-modes) counterparts. However inclusion of 
BLT parameters would shift the masses of KK-modes from their UED values. For keeping our 
analysis simple we stick to the case of two different BLT parameters. The hrst one ry, specihes 
the gauge and Higgs BLTs while the second one, Vf stands for equal fermion and Yukawa BLT 
coefficients. 

Effective interaction for meson decaying to a pair of can be parametrised by a 

4-fermion interaction. The coefficient for this effective interaction is calculable in the frame¬ 
work of SM and in the model of our interest namely the UED model. There are two sets of 
diagrams contributing to this decay process. The Z-penguins contribute dominantly while the 
box diagrams are sub-dominant. Diagrams 3, 4, 6, 7 and 8 in Fig. |2] and Ib and Hb of Fig. [3] 
captures the effects of top-Yukawa couplings, due to which this one loop mediated process is 
being amplihed. 

We have listed the amplitudes from each diagram separately in Appendix A and B. The total 
contribution coming from the penguins are hnite and GIM mechanism plays a crucial role to 
tame the divergences. It should be noted that the hnal results is sum of contributions coming 
from different KK-levels and it also contains the SM contribution (i.e. from 0*^ KK-mode). 
While summing over the KK-levels we restrict ourselves to 5 levels in view of a recent analysis 
relating the Higgs boson mass and cut-off of a UED theory [29]. It has been shown explicitly 
that the limits on the parameters derived from our analysis would change a little if we take 20 
KK-levels in the summation instead of 5. 

There is a one (two)-standard deviation difference between the experimentally measured branch¬ 
ing ratio of Bs(d} with its SM prediction. We have used the experimental data on the 

measured value of branching ratio of i?s(d) —)■ to constrain the parameter space of nmUED 

model. As our calculation would reproduce the results of Br{Bs —)■ /i'''/i“) in the vanishing 
BLKT limits, at the very outset, we have used the above algorithm to set a lower limit on R~^ 
in the framework of UED. In case of UED, present analysis constrained R~^ to greater than 
450 GeV at 95 % G.L. This limit in the framework of UED, is not one of the most stringent. 
However, this is in the same ballpark with the limits those are obtained from the consideration 
of Rb [IB] or p-parameters [52] • The nmUED model has more than one parameters apart from 
the compactihcation radius. The BLT coefficients control the masses of KK-modes and their 
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interactions. So in the nmlJED framework, any lower limit on R~^ would depend on these BLT 
coefficients. As for example, for Rf = Ry = 3, R~^ > 0.7 TeV, while for Rf = 4, Ry = 3, 
R~^ > 0.8 TeV. These lower limits on R~^ are so far the most stringent one in the framework 
of nmlJED. Thus the recent experimental result of Bg meson decay to combined with 

present analysis has excluded the largest region in the nmUED parameter space. Unfortunately 
the bounds on the parameters using the Br{Bd —t is not so competitive. Lower limits 

on R~^ derived from our analysis for negative values of Rf are not so stringent and the have 
been already ruled out from the consideration of electroweak precision data. 

We have also reviewed the effect of electroweak precision constraint on nmUED parameter 
space. Our analysis reveals that the experimental data on the Bg —)■ branching ratio 

would put more severe constraint on the lower limit of R~^ for positive values of Rf than the 
electroweak precision data. 
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Appendices 

A Different contributions to Z^-penguin and self-energy 
diagrams 


The contributions of the Z-penguin diagrams to the functions in Fig. |2] are given 


as follows: 
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where the functions hg and are given by: 
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(A-10) 


The contributions of the self-energy diagrams to the functions Cn are given by the following 
expressions which are generated from Fig. [S] 
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Here A = ^ -|- In dvr — ye, D = 4 — e. 


B Different contributions to box diagrams 


The contributions of the box diagrams to the functions Bn in Fig. [Hare given as follows: 
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The functions U and U are defined as, 
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C Feynman rules for iii nmUED 

The Feynman rules for the different vertices with the assumption that all momenta and helds 
are incoming. 

1) ^ - daiyC, where C is given by: 
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Here the scalar helds S = H,G. 
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5) S^fj2 


92 


\/2Myi( 


n) 


{PrCl + PrCr), where Cl and Cr 


are given by: 
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Here the fermion fields f = u,d, T/, T/, u, e, C^, Ce, S- 
6) W^^Jj 2 ■ ^'y^PLCL, where Cl is given by: 
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The mass parameters ttIx^ are given by: 
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where mj represents the mass of the 0-mode up-type fermion and Cin 
sin(Q!j„) with as dehned earlier. 

And the mass parameters are: 
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where rrij represents the mass of the 0-mode down-type fermion. Here, /" 
integrals given in Eqs. [23] and |23| respectively. 
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= cos(ai„) and Sin = 
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and I 2 are the overlap 
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